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ABSTRACT 

We report multi-wavelength power spectra of diffuse Galactic dust emission from BLAST observa- 
tions at 250, 350, and 500 fim in Galactic Plane fields in Cygnus X and Aquila. These submillimeter 
power spectra statistically quantify the self-similar structure observable over a broad range of scales 
and can be used to assess the cirrus noise which limits the detection of faint point sources. The advent 
of submillimeter surveys with the Herschel Space Observatory makes the wavelength dependence a 
matter of interest. We show that the observed relative amplitudes of the power spectra can be related 
through a spectral energy distribution (SED). Fitting a simple modified black body to this SED, we 
find the dust temperature in Cygnus X to be 19.9 ±1.3 K and in the Aquila region 16.9 ± 0.7 K. 
Our empirical estimates provide important new insight into the substantial cirrus noise that will be 
encountered in forthcoming observations. 

Subject headings: submillimeter — ISM: clouds, cirrus — balloons 



1. INTRODUCTION 

Energy is being injected continually into the interstel- 
lar medium (ISM) through spiral shocks, violent out- 
flows from massive protostars, stellar winds, expanding 
H II regions and supernova explosions. The result is a 
turbulent medium where the dust is well mixed in the 
structures produced. The emission from the relatively 
nearby ISM at high Galactic latitude is known as Galac- 
tic Cirrus. Cirrus-like structure in the brighter emis- 
sion near the Galactic Plane has be en called 'interstel- 
lar froth' (| Waller fc Boulanger1ll994t ): this might be an 
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interesting distinction with a physical basis, but here 
we will simply use the term 'cirrus'. The dynamics of 
the ISM appear to have made the distribution of den- 
sity structures self-similar, with fluctuations in column 
density and surface brightness present on all observable 
scales, though decreasing towards smaller scales. Statis- 
tical description of random fluctuations through struc- 
ture functions is an important method of extracting phys- 
ical properties hidden in diffuse emission as well as for 
providing a quantitative measure to compare with sim- 
ulations. A common statistical tool used to estimate 
the level of cirrus noise is the power spect rum. This is 
valid f or Gaussian random fields. However. I Gautier et al.l 
(1992) fou nd evidence for non-Gaussian ity and recent 
studies by IMiville-Deschenes et al.l (|2QQ7l ) have revealed 
non-vanishing skewness and excess kurtosis in the un- 
derlying brightness fluctuation fields. Nevertheless, for 
estimating the variance the power spectrum is still in- 
dicative. 

Using power spectra, iGautier et al.l (|1992[ ) quantified 
how fluctuations associated with Galactic cirrus are a 
source of confusion noise which limits the detectability 
of point sources. Even with the improved angular res- 
olution of imagers like SPIRE and PACS on the Her- 
schel Space Observatory 16 (hereafter Herschel), this 'cir- 
rus noise' remains important, often dominant. Therefore, 
to assess the noise it is vital to know the statistical prop- 
erties of the interstellar diffuse emission at the relevant 
wavelengths. 

A number of other statistical analyses of emission have 
been carried out on the basis of power sp e ctra for high 
latitud e clo uds, e.g., by [K iss et al. (2003), Ueong etall 
(|2005D . and IMiville-Deschenes etal] (|2007f ). In this pa- 
per we analyse diffuse dust emission in the Galactic 
Plane and present for the first time mult i- wavelength 

16 http://www.esa.int/SPECIALS/Herschel/index.html 
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power spectra in the submillimeter, based on observa- 
tions in Cygnus X (Cyg X) and Aquila with BLAST 
(Pas cale et all 12008( 1 at 250, 350, and 500 //m. We also 
analyse the same regions at 100 and 60 /xm using IRIS, 
IRAS data reprocessed by IMiville-Deschenes fc Lagacha 
(|2005l ). Compared to high latitude studies, analysing dif- 
fuse emission in the Galactic Plane in terms of the struc- 
ture of the ISM resulting from its turbulent properties is 
more challenging because of the long path lengths, high 
column density, high star- formation rate, and contami- 
nation by compact sources. In addition, in the Galactic 
Plane, self-gravity can play an important role in shaping 
diffuse structures at smaller scales. Of particular impor- 
tance are the potentially strong spatial variations of the 
radiation field due to star formation activity. As dis- 
cussed further in § 16. 1\ this could produce variations in 
the dust emission independently of any changes in struc- 
ture of the ISM. 

Our paper is organized as follows. We begin with a 
brief description of BLAST observations (§ [2]) and then 
introduce important aspects of the power spectrum and 
accompanying cirrus noise (§[3j). We analyse IRIS data 
in § HI placing this in the context of earlier studies and 
providing a short-wavelength reference for our submil- 
limeter studies. In § [5] we examine the BLAST data: 
the noise; effect of the point spread function (PSF) or 
beam; removal of compact sources; and the exponents 
and amplitudes of the submillimeter power spectra. We 
estimate the cirrus noise for these maps and compare 
it with the completeness depth of t he BLAST05 point- 
source catalogs at 250 jira. In § 15.61 we also discuss some 
implications for related approved observations with Her- 
schel. We show in § [6] how the observed wavelength de- 
pendence of the amplitude of the power spectrum can 
be understood as a straightforward consequence of the 
spectral energy distribution (SED) of the dust and we fit 
a simple modified black body to estimate the dust tem- 
perature. Our empirical results provide new insight into 
what cirrus noise might be expected in submillimeter ob- 
servations. 

2. OBSERVATIONS 

In 2005, the Balloon-borne Large Aperture S ubmil- 
limeter Telescope (BLAST05: IPascale et alll2008f ) made 
unbiase d surveys in targeted regions of the Galactic 
Plane. iChapin et all (|2008l ) analysed a 4 deg 2 field in 
Vulpecula and here we used data from the two largest 
surveys. Aquila (Rivera-Ingraham et al., in preparation) 
is a 6 deg 2 region observed for 6.1 h, while Cyg X (Roy 
et al., in preparation) covers 10 deg 2 observed over 10.6 
h. Both have good cross-linking from orthogonal scan- 
ning. The maps were made with the SANEPIC algo- 
rithm (jPatanchon et a l. 2008) and were calibr ated using 
the procedure discussed in lTruch et al.l (|2Q08h . 

Although a 2-m telescope, BLAST05 produced maps 
of only 3' resolution (see § !5.2p due to an anomalous PSF, 
corrupted by some uncharac terized combinatio n of mir- 
ror distortion and de- focus (jTruch et al.l 12008). 17 Nev- 
ertheless, maps from the 2005 flight have high signal-to- 
noise and are oversampled with 15" pixels, so that Lucy- 
Richardson (L-R) deconvolution can be used to improve 

17 This problem was fixed for the 2006 flight from Antartica 
(ITruch et al.ll2009h . 



the resolution significantly (Roy et al., in preparation). 
This goal is particularly important for extracting point 
sources (§ 15.3ft but otherwise not essential for the study 
of diffuse emission. We analysed the two largest surveys, 
selecting in Aquila a square sub-field of size 1.83° cen- 
tred on I = 45.85°, b = —0.12° and a similar-sized field in 
Cyg X centred on I = 79.89°, b = 0.47°. The latter does 
not include the brighter star- f orming regions to the east 
containing W75-N and DR21 ([Schneider et al.ll2006T ). 

We also analysed IRIS maps at 100 and 60 fim toward 
these selected regions. These maps have ~ 4' resolution 
on 1.5' pixels (Miville-Deschenes et al. 2002) and we used 
the version in which the s ources have been removed by 
the te chnique described bv IMiville- Deschen es fc Lagachd 
([2006h . 

3. POWER SPECTRUM AND CIRRUS NOISE 

The power spectrum is the Fourier transform of the 
auto-correlation function of the intensity map I(x, y). In 
fc-space, P(k x , k y ) is simply related to an image by 

P{kx->ky) — (I{kx->ky)I (kxiky))i (1) 

where I(k x ,k y ) is the Fourier transform of the image 

and J* its complex conjugate. We used the IDL routine 
FFT to compute the two-dimensional Fourier transforms. 
To ensure a smooth periodic boundar y condition near 
the edges (jMiville-Deschenes et al.ll2002f h the maps were 
apodized by a sine function over a range 10% the width of 
the map. The power spectrum P(k) is obtained by aver- 
aging P(k x ,k y ) over an annulus placed at fc = ^jk\ + fcj*. 

In practice, contributions to the total power spectrum 
come not only from diffuse dust emission, but also point 
sources, the cosmic infrared background (CIB), and the 
noise. When these components are statistically uncor- 
r elated, the total power spect rum can be expressed as 
(jMiville-Deschenes et aTll2007[ ) 

p(fc) = r(fc) [P cirrus (fc) + 

-^source 

(k)+P C m(k)] + N(k), 
(2) 

where T(fc) is the power spectrum of the PSF (the square 
of the modulus of the two-dimensional Fourier transform 
of the PSF), which decays at large fc. For the bright 
Galactic Plane fields targeted here, the contribution from 
the CIB to the power spec trum is negligible. While the 
noise is measurable (§ 15. 5|) . it too makes an insignificant 
c ontribution. 

iGau tier et al. (1992), followed by many other author s 
(jKiss et all I2QQ1L 120031 : IMiville-Deschenes et al.l I2002D . 
have shown that the power spectrum of Galactic cirrus 
follows a power law 

P(k)=P(k )(k/k ) a , (3) 

quantified by an amplitude Po = P(fco) at some 
fiducial fcp and an exponent a that is typically —3 
(Miville-Deschenes et al. 2007). From their analysis, the 
associated 'cirrus noise' for a telescope with mirror di- 
ameter D working at wavelength A can be quantified as 

---"»<'/">"(T^£f (4, 

[P(fc = 0.1 arcmin-^/lO" 3 MJy 2 sr -1 ] " 5 mJ y- 
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Here our assumptions about the measurement strat- 
egy for point sources are the same as made by 
Helou & Beichman (1990) and a dopted bv IKiss et al.l 
(|2001f ). which in the notation of IGautier et al.l (|1992) 
are a measuring aperture d, with 'resolution ratio' r = 
d/(X/D) of 1.6 and a reference annulus with 'separation 
ratio' 2, and also a = —3, close to what we find below. 
The beam-related factor quantifies the effect of a smaller 
beam probing smaller spatial scales where the power in 
the fluctuations is weaker. 

If the sources being measured are extended, as will be 
the case in many Galactic surveys, then the measuring 
aperture needs to be larger. The main consequence of 
this on increasing cr C i rrus is captured by t he factor r 2 5 
for th e range of interest (see also Fig. 3 in IGautier et al.l 
fl992h . 

Working with 100 fim IRAS data on fields of differ- 
ent average surface brightness (2ioo)> an d adopting a 
fiducial scale kp = 0.01 ar c min - 1 at which the ampli- 
tude is Pioo> Ga utier et al.l ([1992) found the trend that 
-Ploo = C(hoo) 3 , where C is a proportionality constant. 
If this is substituted in the abo ve, we recover the formula 



given by Helou fc Beichmanl (|1990l ) and evaluated by 
IKiss et al.l (|2001h estimating the cirrus noise for different 
levels of cirrus bright ness. Note that the cirrus no ise es- 
timate described bv IMiville-Deschenes et al.l (|2QQ7l ) used 
a slightly different definition of the noise, as well as in- 
corporating a trend giving a slightly steeper power law 
with increasing (2ioo)> an d in effect a lower C, which to- 
gether conspire to lower the noise estimate by about a 
factor six for bright cirrus. 

For reasons discussed below, both estimating P100 from 
(iioo) and scaling it to the appropriate wavelength of ob- 
servation are problematical. Therefore, if possible the 
exponent and amplitude of P{k) should be measured di- 
rectly for the field of interest at the relevant wavelength. 
Measuring and normalizing at a scale as close as possible 
to the beam also avoids issues of extrapolation. With 
Herschel in mind, we have chosen k = 0.1 arcmin -1 . 
The normalization of the amplitude in equation (|4]) an- 
ticipates what we find at 250 fim for the two bright fields 
we examined. 

4. IRIS POWER SPECTRA 

Our primary goals are to analyse statistical fluctua- 
tions of diffuse dust emission in the Galactic Plane us- 
ing the power spectrum, and to measure the wavelength 
dependence of the amplitude of this power spectrum. 
To connect with previous work and provide a reference 
at shorter wavelengths, we begin with IRIS. IRIS (like 
IRAS) data come in 'plates' 12.5° on a side. We com- 
puted the power spectrum for the entire plate containing 
the Aquila field (plate 263) but because of a gap in the 
IRAS sky coverage near Cyg X (plate 361) we had to 
settle for a smaller region, 6° on a side in Galactic co- 
ordinates centred on I = 80.69° and b = 0.76°. Binned 
estima tes o f P(k) and their uncertainties were evaluated 
as in § 15.41 

The power spectrum for the Aquila field at 100 /xm is 
shown in Figure [TJ the power spectrum for the Cyg X 
field (not shown) is similar. The noise, N(k), can be es- 
timated by comparing different independe nt IRAS maps 
(HCONs; IMiville-Deschenes et all l2002D and for these 
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Fig. 1. — IRIS 100 /im power spectrum (dot-dash curve) of a 
square 12.5° on a side in Aquila spanning the Galactic Plane; and 
after dividing by T(k) (solid). The error bars, plotted only every 
tenth point for clarity, do not account for the imperfect approxi- 
mation to the beam whose effect is clearly present at large k. We 
fit the power law (dashed line) over the range 0.004 arcmin -1 < k 
< 0.1 arcmin -1 marked by the vertical dashed lines. The exponent 
is -3.05 ± 0.02 and P 100 = 0.63 ± 0.03 MJy 2 sr -1 . 

bright fields is negligibly small. The behavior of P(k) 
is very much like that seen at higher latitudes and with 
BLAST (§ I5.5|) , except that here it does not fall off at the 
highes t k quite as fast as expected f rom the effect of the 
beam ([Miville-Deschenes et al.f2002l ). We have examined 
a number of plates crossing the Galactic Plane as well as 
the adjacent plates at higher latitude and conclude that 
this is an effect seen only in the Plane. Apparently in 
these source-removed maps there are significantly larger 
fluctuations within the scale of the beam than expected, 
and subsequently these are greatly amplified by the beam 
correction. These fluctuations are likely an artifact of the 
source removal technique, which would be maximized in 
the Galactic Plane, where the surface density of sources is 
greatest . For the beam r(fc) we used a Gaussian as spec- 
ified by IMiville-D eschene s et al.l (|2002l ) which is clearly 
just an approximation to the complicated effective beam 
of the images formed from the IRAS timest reams. In 
any case, here we use only the part of the power spec- 
trum at lower k < 0.1 arcmin -1 , which is largely un- 
affected by the beam correction. We also avoid small 
values of k < 0.004 arcmin -1 where the power is signif- 
i cantly affected by the apodiz ation of the image edges 
(jMiville-Deschenes et al.l l2002f ) . Figure [T] confirms em- 
pirically that this is a good choice for the range of k used 
for the power law fit. The non-linear weighted fits were 
carried out with the IDL routine. 

For 100 /im we find an exponent a of —2.71 ± 
0.07 in Cyg X an d -3 .05 ± 0.02 in Aquila. 
IMiville-D eschene s et al.l (|2007| ) find that the steepness 
of the power spectrum increases as a slow function of 
(hoo). For the two fields (I 100 ) is 260 and 120 MJy sr -1 , 
respectively, and the expected exponents from this trend 
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are —3.4 and —3.3, slightly steeper than what we find. 
Note however that they find a scatter of 0.3 about their 
trend. Note also that while we report the formal errors 
from the fit to the power spectrum, systematic errors in 
the derived exponent could be larger, on the order of 0.1. 

We find Pioo = 1.61 ± 0.22 and 0.63 ± 0.03 MJy 2 sr" 1 
for Cyg X and Aquila, respectively. These can be com- 
pa red to values of 4.3 and . 59 MJ y 2 sr _1 from the trend 
in iMiville-Deschenes et al.l (|2007[ h about which there is 
a factor of three scatter at high (Jioo)- Again the agree- 
ment is better for the Aquila field, despite the fact that 
there are significant asymmetries in the image (the bright 
swath of the Galactic Plane and residual striping) which 
make the two-dimensional power spectrum not quite cir- 
cularly symmetric. 

For more direct comparison with the BLAST power 
spectra, we computed the IRIS power spectra for the 
smaller BLAST sub-fields. The exponents of the power 
spectra are —2.97 ± 0.23 and —2.82 ± 0.10, respec- 
tively. These are harder to measure, given the smaller 
dynamic range, but do not appear to have changed 
significantly despite the larger (Jioo) values (430 and 
330 MJy sr _1 ). Expressing amplitudes as Pioo (though 
k = 0.01 arcmin -1 is beyond our range in these small 
sub-fields), we find 3.1 ± 1.3 and 0.75 ± 0.11 MJy 2 sr" 1 , 
respectively, compared to trend values of 29 and 13 
MJy 2 sr _1 . Because of the P 1 / 2 dependence in equa- 
tion (|4]), predictions of cirrus noise depend less strongly 
on any deviations from the trends. Nevertheless, these 
results illustrate the important point that when attempt- 
ing to assess the cirrus noise, one should, if possible, 
measure the exponent and amplitude near the spatial 
frequencies of interest. 

Because the cirrus noise is wavelength dependent, ide- 
ally Po would be measured at each relevant wavelength 
too. The amplitude might in principle be scaled, say from 
Pioo- Not all scaling prescriptions in the literature can 
be valid, however, and we discuss our recommendation 
in §E1 

5. BLAST POWER SPECTRA 

5.1. Noise Power Spectrum 

In the BLAST map-making procedure 
(jPatanchon et al.l l2QQ8h . a variance map V(x,y) is 
produced based on the noise in the time-stream data 
and the map coverage by the bolometers in the arrays. 
For these BLAST surveys, the resulting maps of V are 
quite uniform. A realization of the noise map N(x,y) 
can be constructed for each independent, uncorrelated 
pixel from N(x,y) = ^/V(x, y)v, where v is a Gaussian 
random variable with unit standard deviation and zero 
mean. Given the relative uniformity of the variance map, 
the resulting simulated noise map is close to white on all 
scales and thus the power spectrum of this noise is quite 
flat (see the example for Cyg X at 250 jam in Fig. [2|). 
At large k the total power spectrum decays to this noise 
level, because of the combined effects of the decreasing 
power in the diffuse cirrus emission and the PSF. There 
is also a component of low frequency noise in the map 
arising from the very long time scale 1/ f noise present 
in the time streams. However, in cross-linked maps 
produced using data from multiple scanning directions 
this is greatly reduced by the SANEPIC algorithm (e.g., 




0.01 0.10 1.00 

k [arcmin -1 ] 



Fig. 2. — BLAST power spectra at 250 /im for the Cyg X region 
shown in Fig. 0] The dotted curve, corresponding to the upper 
panel of Fig. 0] shows that a significant amount of power is present 
at intermediate scales due to compact sources. The lower dot-dash 
line is P(k) of the same region after removing sources (correspond- 
ing to Fig. [5}. The oscillatory behavior at large k is due to the 
corrupted BLAST05 PSF (Fig. [3}. The almost horizontal dashed 
line is the power spectrum of the noise map (§ I5.1|) . At very small 
scales where the astronomical signals become correlated within the 
beam, the power spectra meet the noise level. The solid line is 
the power spectrum P c irrus(^) 5 after subtracting the noise and di- 
viding by r(/c) to remove the effect of the beam (as in Fig. [T] 
error bars at high k do not include the imperfect knowledge of the 
beam and noise). The vertical dashed lines indicate the restricted 
range 0.025 arcmin -1 < k < 0.2 arcmin -1 used for fitting the 
power law. The exponent is —2.60 ±0.07 and P(0.1 arcmin -1 ) = 
(1.60 ±0.05) x 10 -3 MJy 2 sr -1 . 

Fig. 10 in iPatanchon et aD l2QQ8h and is not important 
here because the cirrus signal at small k is so large. 

5.2. Effect of the Beam 

The BLASTQ5 PSF in telescope coordinates is shown 
in iTruch et al.l ([2008). A synthetic beam can be made 
for a particular map, taking into ac count the scan an- 
gles and coverage (|Chapin et al.ll2008l ) . Because this syn- 
thetic beam is not known out to the full size of the map, 
we have not derived its power spectrum T{k) directly. In- 
stead, we convolved the noise map with the PSF directly, 
found the power spectrum of that map, and divided by 
the power spectrum of the noise map. Figure [3] shows 
r(jfe) for the BLAST05 PSF at 250 /im for the Cyg X 
field. This has a number of features characterizing scales 
seen in the corrupted PSF. Clearly, the observed P(k) 
will be seriously suppressed at high k and so we will use 
data only for k < 0.2 arcmin -1 . At these scales, where 
the beam correction is not too large and so both more 
reliable and of less import, T(k) can be described by a 
Gaussian with <Jr = 0.08 arcmin -1 . Using the Fourier 
relation 2y/27rar cr\) = 1, this corresponds to a Gaussian 
beam of <r b = 1.4', or FWHM = 3.3', slightly smaller 
than for IRIS. T his is close to the 3 .1 f full width at half- 
power found bv ITruch etHI (|2008f ). The corruption of 
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Fig. 3. T(k) for BLAST05 in Cyg X at 250 /an (dashed curve). 
The dot-dash line shows how this can be approximated by the 
power spectrum of a Gaussian PSF of FWHM 3.3'. 

the BLAST05 PSF is such that the initial falloff in T(k) 
is rather similar for the three wavelengths, and for both 
fields. 
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80.0 79.0 
Galactic Longitude 

Fig. 4. — BLAST05 map at 250 of the selected region of 
Cyg X used for studying diffuse emission. The map in the upper 
panel includes point sources which, due to the corrupted BLAST05 
beam, are compact structures of size ~ 3.3 X . In the lower panel the 
sources have been removed. 



5.3. Source Removal 

In the Galactic Plane, the power spectrum is seriously 
contaminated by compact sources. To remove them, we 
first deconvolve the BLAST maps using the L-R method 
which, importantly, conserves flux. Compact sources are 
of size ~ 1.5' in these maps and more easily identified. 
They are fairly well described by Gaussians. We fit Gaus- 
sians (multiple if crowded) to obtain flux densities, posi- 
tions, and FWHM (major axes and position angle) of the 
compact sources. We then convolve these Gaussians with 
the synthetic beam and subtract them from the original 
maps. For the Cyg X sub-field analysed, the upper and 
lower panels of Figure [4] show BLAST maps at 250 /xm 
before and after source removal, respectively. Faint resid- 
uals appear near some of the brightest sources, because 
of multiscale structure in the ISM in which they are em- 
bedded. 

5.4. Two-dimensional Power Spectrum 

Figure [5] shows the two-dimensional power spectrum 
for the Cyg X region at 250 /mi. Because a ~ — 3 pro- 
duces a large intrinsic dynamic range in the power spec- 
trum, we have multiplied it by (k/ko) 3 , with ko = 0-1 
arcmin -1 . The dark rings are produced by the same fea- 
tures giving the dips in r(fc). Noise is being amplified at 
large k because of the k 3 multiplier. 

The two-dimensional power spectrum is circularly sym- 
metric, justifying the annular averaging to find the mean 
P(k). Most of the information relevant to the diffuse 
emission is contained in the plateau in the central region 
of Figure El and the average there in this representation 
is close to P(k = 0.1 arcmin -1 ). 




0.00050 



0.00038 



0.00025cx 



0.00013 



0.0 

-1 1 

k x [arcmin -1 ] 

Fig. 5. — Two dimensional power spectrum multiplied by (k/ko) 3 
(with ko = 0.1 arcmin -1 ), for the Cyg X map with compact 
sources removed, as shown in the lower panel of Fig. 0] Dark 
rings are an imprint from the BLAST05 PSF (Fig. |3}. Scaling 
by k 3 has amplified noise at large k and highlighted the bright 
central plateau, which contains the information most relevant for 
this study. 
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The error in P(k) for each annulus is the standard 
error of the mean. While removing the effect of beam 
from the power spectrum, we have only scaled this error 
by the inverse beam, not accounting for the fact that its 
shape is not perfectly known at high k. However, this is 
of no consequence because we do not fit the data at large 
k where the beam correction is significant. 

5.5. Exponent and Amplitude 

Figure [2] shows P(k) for 250 jam for the Cyg X region, 
for the original map and with sources removed. The 
effect of the beam is dramatic at high fc, showing that it 
is important to have a good estimate of the PSF. There 
is also clearly excess power due to sources and so it is 
important to remove them carefully. 

At small /c, the power spectrum is dominated by the 
dust structures; it decays toward higher /c, as seen pre- 
viously in IRIS data at shorter wavelengths. The os- 
cillatory effect of r(fc) can be seen clearly at k > 
0.2 arcmin -1 and at higher k the power spectrum meets 
the independently-predicted noise level shown. 

-Fcirrus(fc)? obtained by dividing the source-removed, 
noise-subtracted power spectrum by T(k), is also shown 
in Figure O It appears to have a power-law form over 
the limited dynamic range in k available. We fit only 
over the range 0.025 arcmin -1 < k < 0.2 arcmin -1 , re- 
stricted for the reasons discussed in § HI In particular, 
at large k the effect of correcting for the PSF is very 
large. The synthetic PSF is an average for the entire 
observed field and so is imperfectly modeled to the pre- 
cision that would be required for precise compensation 
at high k. Convincing evidence for any deviation from a 
power law for the higher k range will have to await the 
higher resolution observations anticipated with Herschel 
which will probe to beyond k = 1. At smaller k the ob- 
served power spectrum is affected both by apodization 
and by the effective filtering of the map-making proce- 
dure, such that the largest scales are not recovered. For 
the relatively small size of these maps and the charac- 
teristics of these BLAST0 5 observations these e ffects set 
in at about the same k (jPatanchon et al.l 12008). In ear- 
lier work with IRAS, the amplitude Po is often cited for 
fco = 0.01 arcmin -1 . This fiducial value is below the 
range of our observations. Using fco = 0-1 arcmin -1 
would seem preferable and probably more relevant to 
small maps at higher resolution to be made with Her- 
schel. 

Figure [6] shows (k / k Q f P ciTTUS (k) for all three BLAST 
bands (250, 350, and 500 /im) for Cyg X. We find 
P(0.1 arcmin -1 ) = (1.60 ± 0.05) x 10 -3 MJy 2 sr -1 at 
250 /im. Po increases with decreasing wavelength for 
the BLAST bands and for the IRIS bands (§ BJ remains 
about the same at 100 /im and then decreases at 60 /im. 
We obtain an exponent a equal to — 2.60±0.07 at 250 jam. 
This can be compared with the value found above at 
100 /im, —2.97 ± 0.23. The exponents are quite simi- 
lar for all b and s, except for 60 jam which is discussed 
further in § 16.11 If the power spectrum were less steep 
moving into the submillimeter, this would provide impor- 
tant evidence for a change in the statistical properties of 
the emission and the underlying mass column density 
and temperature distribution (§ 16. ip . Definitive searches 
for such wavelength dependence should benefit from the 




10" 5 I t t t 

0.01 0.10 
k [arcmin -1 ] 

Fig. 6. — (k/k ) 3 P c i rrus (k) for the three BLAST bands and two 
IRIS bands, for the Cyg X field with ko = 0.1 arcmin -1 . The 
exponents are similar for all bands. The amplitude increases with 
decreasing wavelength for the BLAST bands and IRIS 100 /im band 
(solid), and then decreases at 60 fim (dotted). Lines represent a 
power law with common exponent —2.65 fit to the data plotted 
with open circles. Data at higher k are not used because of the 
large and uncertain beam correction and eventually the dominance 
of noise (see Figs. [3] and [2j. 




0.01 0.10 
k [arcmin -1 ] 



Fig. 7. Like Fig. [6] but for the Aquila region. The amplitudes 
at 100 and 60 /im are relatively lower, which indicates that the 
cirrus in this region is cooler. The common exponent used for the 
fits was —2.91. 

higher spatial dynamic range anticipated in observations 
with Herschel. 

Figure [71 presents the results from our analysis of the 
Aquila field. Overall the same behavior is seen, though 
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the exponent appears to be slightly steeper (— 3.14±0.10) 
at 250 jum, P at 250 jam is larger, (3.03 ± 0.15) x 10~ 3 
MJy 2 sr _1 , and the relative amplitudes at the IRIS bands 
are much lower. 

5.6. Catalog Depth 

Equation (j4j) can be used to predict the depth reached 
in a compact source catalog. To apply this to the 
BLAST05 Cyg X survey, (where the compact sources are 
of apparent size ~ 1.5' in the L-R maps in which they are 
detected and measured), the appropriate r is ~ 3.6. With 
D = 1.9 m and P(0.1 arcmin -1 ) = (1.60 ± 0.05) x 10~ 3 
MJy 2 sr _1 for the somewhat dimmer part of the survey 
field examined here, cr C i rrus = 4.2 Jy at 250 jam. This 
compares well to the empirically-estimated 3-cr detection 
threshold of ~ 15 Jy (Roy et al., in preparation). The ef- 
fective noise near bright sources can be somewhat higher 
because of contaminating artifacts produced in the L-R 
deconvolution. 

For the large survey of Vela carried out by BLAST in 
2006, where the beam of the 1.8 m mirror was diffrac- 
tion limited (Netterfi eld et all 120091 ). we have not mea- 
sured Po directly, because the region was not covered 
by orthogonal scanning and so the cirrus structure is 
not well constrained over all s cales in the cross-scan di- 
rection (|Patanchon et all |2008[ ). Nevertheless, follow- 
ing the scaling suggested in § [6j we can estimate Po 
from (iioo) ~ 100 MJy sr _1 and T d ~ 18 K, giving 
P(0.1 arcmin -1 ) = 10 -3 MJy 2 sr -1 at 250 /mi, from 
which cr C i rms ~ 0.6 Jy. However, the sources were actu- 
ally extended, with typical apparent sizes of 1', so that 
r ~ 2.4, and then cr c i rrus ~ 1.6 Jy. For comparison, the 
catalog d epth judged from simu lations of completeness 
was 6 Jy ([Netterfield et al.ll2009f ). We conclude that the 
depth is dominated by the influence of the cirrus noise. 

With Herschel and assuming the source sizes are com- 
patible with r = 1.6, cr C i rms in the Cyg X region should 
be closer to 100 mJy at 250 jam. While a distinct im- 
provement, this is nevertheless substantial compared to 
the instrument noise (r.m.s. ~ 5 mJy) predicted using 
HSPOT 18 for the planned parallel PACS-SPIRE map 
scanning strategy in the guaranteed-time key project 
HOBYS 19 (Herschel imaging survey of OB Young Stellar 
objects). We expect substantial regional variations. For 
example, in Aquila, simply from the larger Po the depth 
achieved should be 1.4 times worse. And this is not even 
the brightest part of the Galactic Plane to be surveyed in 
the open-time key project Hi-GAL 20 (Herschel infrared 
Galactic plane survey). 

Extraction of sources with wavelength-dependent and 
multi-scale structure will be especially problematical and 
it will be necessary both to measure the cirrus properties 
and to simulate its effects. For SPIRE maps, in partic- 
ular those made with the parallel PACS-SPIRE mode, a 
corollary is that the high signal-to-noise ratio resulting 
from the redundant coverage is well suited to studying 
the statistical properties of the bright cirrus. Of particu- 
lar interest is how the power spectrum varies with wave- 
length, which can be used to determine physical proper- 

18 http://www.ipac.caltech.edu/Herschel/hspot.shtml 

19 http://starformation-herschel.iap.fr/hobys/ 

20 https:/ /hi-gal. ifsi-roma.inaf.it/higal 



ties of the diffuse dust. 

6. WAVELENGTH DEPENDENCE 

6.1. Mass Column Density Distribution 

BLAST and IRAS maps record surface brightness / 
which, in turn, depends on the dust mass column density 
M d : 

h = M d K u B u (T d ), (5) 

where k v is the dust emissivity, B v is the Planck function 
for dust temperature T d and it is understood that the 
right-hand side is summed over various dust components 
which might have different n and T^. 

The underlying spatial property being probed is M d 
which is in turn the projection of the three-dimensional 
density distribution. M d is modulated as a function of 
frequency by the emissivity and temperature, and can 
also have modulations caused by spatial differences in 
these properties. The ISM in the Galactic Plane is far 
from homogeneous, containing diffuse atomic and ionized 
regions and shielded molecular regions with possibly dif- 
ferent mixes of dust compositions and sizes, and along 
any line of sight there will be a range of interstellar ra- 
diation field intensities, all of which could affect T d . The 
surface brightness I v is particularly sensitive to T d for 
frequencies near or above the peak in the SED, i.e., for 
100 jam and shorter for typical interstellar diffuse dust 
temperatures. Therefore, it is somewhat surprising that 
the power spectra exponents that we obtained are so sim- 
ilar. The potential for wavelength dependence should be 
particularly high at 60 jam where non-equilibrium emis- 
sion from very small grains (VSGs) starts to become 
important (jDesert et all 119901 : lLi~fc Drainell2001f ) and is 
even more directly responsive to the ultraviolet radiation 
field. The Cyg X field analysed definitely has different 
spatial structures appearing in the images at the shorter 
wavelengths and there is some evidence in Figure [6] for an 
effect on the exponent of the power spectrum. Further- 
more, the exponent for a particular region must depend, 
through the structure of M d , on the environments probed 
along the line of sight and the energy injection at large 
scales that is responsible for the apparent turbulence in 
the ISM. 

In the simple case of homogeneous correlation in the 
ISM, I v everywhere in the map would scale with fre- 
quency according to the simple relative spectral en- 
ergy distribution S v of the emitting dust. From equa- 
tion ([5]), S„ oc KjjBjj^^). Therefore, Po, measuring 
the same structure in M d , would scale simply as S% 
(see eq. pQ) while cr c i rms would scale as S v (see eq. [3]). 
T his scaling is wha t is a ssumed by the prescriptions 
of Lagac he fc Pugetl (|200Q[ ) and iMiville-Deschenes et al. | 
(120071) , for example, and appreciated bv Ueong et al.l 
(2005;[). However, on making the substitution Pioo oc 
(Jinn) 3 , iHelou fc Beichmanl (|1990l ) end up ass uming that 
P sca les as S 3 and thus cr C i rrus scales as S 1 ' 5 . iKiss et all 
(2001) find some empirical support for the latter over the 
ISO ISOPHOT wavelength range. 

Note that the wavelength dependence of the signal- 
to-noise ratio of a point source will vary by the factor 
(D/A) 2,5 from equation (|2j), and also the ratio of the 
SED of the source relative to that of the confusing cir- 
rus. At submillimeter wavelengths, the situation is more 
favorable for detecting cold sources in warm cirrus than 
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vice versa. Therefore, the completeness depth of a sur- 
vey as a function of wavelength will depend on these 
factors. Sour ce size and structure can also change with 
wavelength (jNetterfield et al.1l2QQ9l ) but we will not dwell 
on this quantitatively here. 

6.2. Characteristic Temperature 

The similarity in the exponents of the power spectra 
as a function of wavelength is no guarantee of underly- 
ing homogeneous conditions. Nevertheless, we can as- 
sume that there is some characteristic SED for each re- 
gion, and use the relative frequency dependence of P^ 2 
to recover it. To be consistent with this picture , we 
evaluated P(0.1 arcmin -1 ) using a common exponent, 
which is simply the average of the exponents found at 
each wavelength (excluding 60 jam). This relative SED 
is shown in Figure [8] for each region. 

Fitting this relative SED with a modified black body 
(eq. [5]) yields a characteristic dust temperature. We 
have assumed k v oc v" with (3 = 2 (Hildebra53Ii2l3). 
We did not include the 60 jam data, because the excess 
there, presumed to be from non-equilibrium emission, 
would be inconsistent with a simple single-temperature 
modified black body. For Cyg X, the increase at long 
wavelengths with respect to the fit is probably the result 
of the range of temperatures known to be present in dif- 
ferent parts of this region (Roy et al., in preparation). 
Therefore, even in this case we kept (3 = 2, although 
empirically (3 = 1.5 would produce a better fit. 

For Cyg X, the best fit temperature is 19.9 ±1.3 
K, which appears to be slightly warmer than the local 
high latitude diffuse m edium dust temperature of 17.5 K 
(Boulang er et al.lll996l ). Finding warmer dust in Cyg X 
is not particularly surprising, because this region, which 
lies along the local spiral feature (I = 80°), is one of 
active star formation, with an OB association and nu- 
merous HII regions. In the 21-cm radio continuum this 
region includes areas of strong diffuse thermal emission. 
To account for the higher equilibrium temperature the 
ambient radiation field absorbed in Cyg X would have to 
be about (19.9/17.5) 6 = 2.1 times higher than locally. 

By contrast, for Aquila the characteristic temperature 
derived is somewhat lower, about 16.9 ±0.7 K. This line 
of sight passes through the inner Galaxy (I = 45°), with 
significant molecular clouds but much less star formation 
(Rivera-Ingraham et al., in preparation). 

As we have noted, under uniform conditions, maps at 
different wavelengths will have to be highly correlated in 
their spatial structures, and so they will have the same 
power spectrum scaled by S 2 . But the fact that the 

1 /2 

relative SED derived from P 1 appears reasonable is no 
guarantee that the maps are simply scaled versions of one 
another (since phases could be different). However, this 
can be checked directly. Inspection of the BLAST maps 
shows that they are remarkably similar in all three bands. 
This is quantified through the correlations with respect 
to the 250 /am map shown in Figure [9] for Cyg X. The 
linear fit s shown are the o rdinary least squares bisector 
solution (jlsobe et al.lll990l ) obtained with the IDL rou- 
tine SIXLIN. Note that the correlation of BLAST 250 /am 
with IRIS 100 jam is not as high as it is between BLAST 
bands. This is expected, because of the above-mentioned 
greater sensitivity of 100 jam emission to Differences 



with respect to the submillimeter images become more 
apparent at 60 jam, where in addition the power spectrum 
is found to be somewhat steeper (Fig. [6]). Figure [T0l shows 
the even better correlations for the Aquila region. A clue 
to understanding the good correlation is, of course, the 
very detailed similarity of the power spectra at the dif- 
ferent wavelengths. 

As a consistency check, the slopes of the correlations 
can be used to construct another relative SED. As seen in 
Figures [9] and [lOl these look very similar to the respective 
SEDs based on the power spectra (Fig. [8]) and the best- 
fit temperatures 18.6 ± 0.5 K and 17.3 ± 0.3 K are very 
similar to those found above. 

We therefore have some evidence that the appropriate 
scaling of a c [ rrus is as 5, not 5 1 ' 5 . However, this is not 
immediately generally useful, unless there is some inde- 
pendent evidence about the frequency dependence of S 
(or the effective and /3), since this can clearly vary 
significantly from region to region. The best approach is 
to measure Po directly at each wavelength. 

7. CONCLUSION 

In these fields in the Galactic Plane, the expo- 
nent of the power spectrum of the 100 jam IRIS 
ma p is close to —3, w i thin the dispersion seen 
by IMiville-Deschenes et al.l J2007) despite the average 
brightness (Iioo) being beyond the range studied by 
these authors. On the other hand, the amplitudes of 
the 100 /am power spectra estimated for these bright 
~ 2°x 2° fields are significantly below what would 
be extrapolated from the tre nd with (iioo) found by 
IMiville-D eschene s et al.l (|2007[ ). Therefore, particularly 
in such bright star forming regions, it is recommended 
that the power spectrum be computed directly. The 
power spectra derived from the BLAST observations are 
also well fit by power laws, with similar exponents. The 
frequency dependence of the amplitude of the power 
spectrum can be described by the square of an SED which 
is a simple modified black body function with a reason- 
able characteristic temperature. This is confirmed by 
direct correlations between the maps at different wave- 
lengths. However, this characteristic temperature does 
appear to change in different Galactic environments, and 
unless its value is known independently the power spec- 
tra and/or map correlations need to be evaluated for all 
wavelengths. 

Cirrus noise will be important in many planned multi- 
wavelength Galactic Plane, high latitude, and extra- 
galactic surveys carried out with Herschel. Our results 
provide important empirical support for a proposed pre- 
scription for the wavelength dependence of the cirrus 
noise, which incorporates a factor which varies directly 
as the SED. In practice the noise will be best evaluated, 

most free of assumptions, simply by measuring P^ 2 at 
each wavelength and using equation (|4]). 
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Fig. 8. — Relative SED obtained from the square root of the amplitudes of the power spectra fit in Figs. [6]and[3 Left Cyg X; right, 
Aquila. Solid curve is the best-fit modified blackbody with (3 = 2. The 60 /im power spectrum values (open circles) were not used for 
the fits (see text). The best-fit temperatures of Cyg X and Aquila are 19.9 ± 1.3 and 16.9 ± 0.7 K, respectively. For ease of comparing the 
two fields and results in Figs. [9] and [10] data and curves have been displayed normalized to the value of the best-fit modified blackbody at 
100 /im. 

(CSBF) staff for their outstanding work. 
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Fig. 9. — Pixel-pixel correlation of 350, 500, and 100 fim maps with respect to the 250 fim map, with the best-fit correlation line plotted. 
The lower right panel shows a relative SED obtained from the slopes of these correlations. With (3 = 2, the best-fit temperature is 
18.6 zb 0.5 K. The normalization here is the same as in Fig. [8] The value at 250 /mi, indicated by the triangle in this normalization, is not 
used explicitly for the fit. The open circle (from the 60 /im slope) was also not used. 
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Fig. 10. Like Fig. [9] but for the Aquila region. The slope of the 100 fim correlation is much shallower than in the Cyg X field, and the 
best-fit temperature is 17.3 ± 0.3 K. 



